Percolation, depinning, and avalanches in capillary condensation of gases in 

disordered porous solids 
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We propose a comprehensive theoretical description of hysteresis in capillary condensation of gases 
in mesoporous disordered materials. Applying mean-field density functional theory to a coarse- 
grained lattice-gas model, we show that the morphology of the hysteresis loops is influenced by 
out-of-equilibrium transitions that are different on filling and on draining. In particular, desorption 
may be associated to a depinning process and be percolation-like without explicit pore-blocking 
effects. 
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Hysteresis in the capillary condensation of gases ad- 
sorbed in disordered mesoporous materials has been 
widely studied, but its theoretical interpretation remains 
controversial. Experimentally, it is mainly observed in 
sorption isotherms measuring the amount of fluid present 
in the solid as the pressure of the ambient vapor (or the 
chemical potential n) is gradually increased and then de- 
creased. Experiments have been performed in various 
adsorbents, from low-porosity glasses to very light silica 
aerogels, and the main observations can be summarized 
as follows: i) the hysteresis loop shrinks as the tempera- 
ture increases and vanishes above a certain temperature 
Th|l|, |^, [|; ii) draining occurs over a narrow range of 
pressures, which usually results in an asymmetric shape 
of the loop; iii) the shape also depends on the geomet- 
ric and energetic properties of the solid and may change 
from smooth to rectangular as the porosity increases J4); 
iv) when the hysteresis is smooth, there is a whole fam- 
ily of scanning curves obtained by performing incomplete 
filling-draining cycles jl], 0, and in the early stages of the 
drainage, the emptied regions exhibit long-range fractal 
correlations, reminiscent of invasion percolation j6|. The 
conjunction of these observations cannot be explained by 
the various mechanisms of hysteresis that have been pro- 
posed in the literature, either modelling the solid as a 
collection of independent pores Q or focusing on kinetic 
network effects (e.g., pore-blocking) 0. In an effort to 
provide a single unifying theoretical framework, we have 
recently presented a new description of capillary conden- 
sation in disordered materials that allows for the study of 
both equilibrium and out-of-equilibrium phenomena^). 
Applying mean- field density functional theory (i.e., local 
mean-field theory or MFT) to a lattice-gas model that 
incorporates at a coarse-grained level the geometric and 
energetic disorder of an interconnected porous structure, 
we have shown that the experimentally observed hys- 
teretic behavior is related to the appearance of a complex 
free-energy landscape with a large number of metastable 
states; the main result is that hysteresis can occur both 
with and without a true equilibrium phase transition. 

In this letter, we use the same theoretical framework 




FIG. 1: Schematic out-of-equilibrium phase diagram of the 
model showing the different regimes discussed in the text. R 
parametrizes the strength of the solid-induced perturbation 
(the units are arbitrary). 



to show that the morphology of the hysteresis loop is af- 
fected by out-of-equilibrium phase transitions (which dif- 
fer on adsorption and desorption) and that the presence 
of the interface between the gas reservoir and the porous 
solid has a dramatic influence on the hysteretic behavior. 
A major finding of our study is that desorption can in- 
volve a percolation process without invoking additional 
kinetic mechanisms such as pore blocking. To understand 
the changes of morphology of the sorption isotherms and 
the out-of-equilibrium collective phenomena, we borrow 
from studies on avalanches and criticality in low-T fer- 
romagnetic materials (i.e., Barkhausen noise) || [lo| 
and on driven interfaces and depinning in systems with 
quenched disorder Jl^, Whereas most of these studies 
are at T=0, MFT allows us to consider finite tempera- 
tures, which are relevant for capillary condensation, and 
to relate such phenomena as avalanches and interface de- 
pinning to the properties of the underlying free-energy 
surface. 

The hysteresis behavior seen experimentally can be 
rationalized by considering the phase diagram plotted 



in Fig. 1, that schematically summarizes the present 
study. More details will be given in a forthcoming 
publication jll). R parametrizes the effective strength 
of the solid-induced perturbation and may depend in a 
complicated fashion on the porosity and on the solid- 
fluid interaction. Within MFT, there are lines of out-of- 
cquilibrium critical points separating a region (1) (low 
T, low R) in which both the adsorption and desorption 
isotherms have a jump that represents a first-order-like 
transition (an infinite avalanche) , an intermediate region 
(2) in which the adsorption is smooth and the desorption 
discontinuous, and a region (3) (high T, high R) in which 
the two branches are smooth and the desorption is essen- 
tially percolation-like. These collective phenomena take 
place below Th at which hysteresis first appears. (The 
equilibrium transition discussed previously^) is not con- 
sidered here.) These sharp out-of-equilibrium transitions 
arc well defined in MFT in which thermally activated 
processes are neglected. In experiments, these phenom- 
ena are observable when the experimental time scale is 
smaller than the equilibration time. The fact that the 
measured hysteresis loops are quite reproducible jl| ^, || 
and closely resemble those predicted by MFT strongly 
supports the existence of such a separation of time scales. 
Note also that avalanche events associated to capillary 
condensation are indeed observed^ ^5| . 

As in previous studies M, we use a disordered lattice- 
gas model whose Hamiltonian is 
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where Tj = 0, 1 and (1 — rji) = 0, 1 denote the fluid and 
matrix occupation variables, respectively, and < ij > 
stands for nearest-neighbor pairs of sites. We consider 
the simplest matrix microstructure in which matrix "par- 
ticles" are distributed randomly over a 3-dimensional 
lattice of N sites. The influence of the matrix is then 
specified by two parameters, the average matrix density 
p m that controls the porosity, and the interaction ratio 
y = w m f /w f / that controls the wetting properties of the 
solid/fluid interface. The results are qualitatively similar 
for more realistic matrix microstructures|I4 . The lattice 



is periodically replicated in all directions. To study sur- 
face effects we consider two different setups, by letting or 
not a slab of width AL X free of matrix sites: this defines 
a region (a "reservoir") that only contains bulk vapor 
(at least for the thermodynamic conditions under study) 
with two planar interfaces in the x-direction (see Fig. 3). 

For a given realization {7/;} of the random matrix, the 
fluid density profile {pi} is obtained from minimization of 
the mean-field grand-potential functional^, which yields 
a set of N coupled equations, 

-p(p+w ff Ej/ilPj+S/C 1 -^)])]-! 



CL 0.4 - 



(b) 




-5.5 -5 -4.5 -4 -6 



-5.5 

n* 




mil 



(2) 



FIG. 2: Hysteresis loops in absence or presence of an external 
interface (p* = p/vjff): (a) y — 0.8, (b) y = 1.2, (c) y = 1.5. 
Note in Figs, a and b the artificial extension of the desorption 
branch in the absence of interface. 



where /3 = l/(ksT) and the sum is over the nearest- 
neighbors of site i. These nonlinear equations are solved 
by a simple iterative method, increasing or decreasing 
progressively p to obtain the sorption isotherms (i.e., 
pf(p,T) — (I/AT) 'J2 i pi) or performing more complicated 
trajectories to obtain the scanning curves. The results 
are illustrated here for a bec lattice, either for a single 
sample of linear size L = 96 (N = 2L 3 ) with AL X = 12, 
or, when the slope of the sorption isotherms is steep and 
exhibit strong variations with L, through a finite-size 
scaling analysis after averaging over several hundreds of 
matrix realizations. 

The shape of the hysteresis loop depends on T, p m and 
y. Here, for the sake of clarity, we only discuss the vari- 
ation with the interaction ratio y while keeping T and 
p m constant (then R = y — 1/2 controls the strength 
of the local random fields in the system^)). Typical re- 
sults obtained for T* — k B T/wff — 0.8 and p rn — 0.25 
are shown in Fig. 2. One first notices that introducing 
an external interface has a major influence on the des- 
orption branch, the adsorption branch being itself unaf- 
fected (see alsojl6]]). This influence is strongest at small 
y and vanishes for y = 1.5: the presence of the surface 
reduces the range of stability of the metastable adsorbed 
liquid (the flat, upper part of the desorption curve in 
Figs. 2a,b), thereby narrowing the hysteresis loop. One 
can also see that the discontinuous jumps present in the 
sorption isotherms for y = 0.8 (region (1) of Fig. 1) are 
replaced by smooth variations for y — 1.2 and 1.5 (re- 
gion (3) of Fig. 1) and that the hysteresis loops become 
asymmetric, as observed in experiments. 

Why is the desorption so sensitive to the presence of 
the external surface? For pores of simple geometry, one 
knows that draining is facilitated by the contact with the 
gas reservoir because of the presence of a meniscus that, 




FIG. 3: Cross-section of a 3-D system of linear size L = 96 
during desorption right at the depinning transition: (a) y = 
0.8, p* = -4.705, (b) y = 1.2, p* = -5.44. The initial inter- 
faces with the bulk vapor are located at x = 12 and x = 96. 
The black regions represent the drained sites {pi < 0.3) and 
the grey dots the matrix sites. 
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FIG. 4: Finite-size-scaling collapse of Apf on desorption for 
y = 1.2 with (j,*, = -5.446, v = 0.89, Df = 2.84 (AL* = L/8). 



at a given /x, can traverse the whole pore at no free-energy 
cost. However, this picture cannot be simply transposed 
to disordered porous materials: this is illustrated by the 
large-?/ case for which the interface has no effect at all. 
The key-concept for understanding the results shown in 
Fig. 2 is that of pinning and depinning of a driven inter- 
face in the presence of quenched disorder. Consider the 
desorption process when slowly decreasing p from p S at, 
the chemical potential at bulk liquid-gas coexistence. At 
first, the interface between the external vapor and the ad- 
sorbed fluid is pinned by the disorder, i.e., it only enters 
superficially in the material and this latter stays macro- 
scopically filled with liquid. As one further reduces p, 
several situations are encountered. For a large enough 
interaction ratio y, bubbles of gas may appear in the bulk 
of the porous medium, and the mass adsorbed decreases 
continuously as these bubbles grow, coalesce, and eventu- 
ally extend over the whole pore space; this phenomenon 
is insensitive to the presence of the surface, as illustrated 
by Fig. 2c. For smaller values of y, no bubbles of gas 
appear in the bulk of the material. This can be inferred 
from the fact that the same system without external sur- 
face does not drain until much smaller values of p are 
reached (Figs. 2a,b). The draining mechanism is thus 
the depinning of the interface with the bulk vapor: for 
a certain chemical potential p c , the driving force on the 
interface is large enough to erase the free-energy barri- 
ers, and the interface sweeps through the whole system, 
which results in the rapid drop seen in the desorption 
isotherm. 

Two different cases can occur as p approaches p c from 
above. For a small interaction ratio (Fig. 2a), the sur- 
face of the growing gas region inside the solid keeps a 
well-defined orientation at long length scales (that of 
the initial interface with the bulk vapor), and, exactly 
at p = p c , the desorption proceeds discontinuously, as 
shown in Fig. 2a|T^]. For intermediate values of y (Fig. 
2b), the growing interface has no well-defined orienta- 



tion; it is rather isotropic and self-similar, and the des- 
orption is rapid (at least in the early stage of draining), 
but continuous, as in invasion percolation. The former 
case corresponds to a compact growth of the gas region 
with a self-affine interface and the latter to a self-similar 
growth, characteristic of percolation: this is illustrated in 
Fig. 3 showing the drained regions where pi < 0.3 (since 
for a large majority of sites, pi is found either lower than 
0.3 or larger than 0.8, drained sites can be conveniently 
defined by this threshold) 14 ). It should be stressed that 
the percolation-dominated draining process which is ob- 
served for intermediate values of y does not require the 
introduction of kinetic constraints such as pore blocking. 
The experimental observations for hexane in VycorQ can 
thus be interpreted without invoking this latter mecha- 
nism. 

One can extract the exponents controlling the critical 
behavior as p — > p c by performing a finite-size scaling 
analysis of Apf(p) — pf(p sa t) — Pf{p)- This is shown in 
Fig. 4 for the percolation regime: there, the value of v 
is compatible with that of 3D percolation, whereas the 
value of the fractal dimension Df is somewhat larger, for 
reasons as yet unclear. A similar analysis can be car- 
ried out in the self-affine regime, and more details will be 
given elsewhere [[HJ. At a critical value of y somewhere 
between 0.8 and 1.2, there is a morphology transition 
from self-affine to self-similar growth, which corresponds 
to a transition from a discontinuous desorption to a con- 
tinuous, percolation-like desorption. This is similar to 
what is observed for domain growth in random magnets 
or in fluid invasion ^M, y playing the role of the param- 
eter controlling the wetting properties of the invading 
liquid. 

Since the filling process is not affected by the presence 
of the interface with the gas reservoir, the jump in the 
adsorption isotherm observed in Fig. 2a has not the same 
origin as that in desorption (note that we had to perform 
a finite-size scaling study to conclude on the existence of 
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FIG. 5: Metastable states obtained by performing desorption 
scanning trajectories in the absence of an external interface 
for y = 1.2. Compare with Fig. 2b. 



a jump in the thermodynamic limit jl4|| ) ■ It corresponds 
to the appearance of a macroscopic, connected network 
of liquid regions in the bulk of the porous medium. In 
MFT, this is again a bona fide phase transition, albeit 
out-of-equilibrium: it is similar to the infinite avalanche 
seen in the Barkhausen noise of low-T ferromagnetic ma- 
terials and discussed in Refs.|], [l0|, O]. One then expects 
that the adsorption isotherm should change from discon- 
tinuous to continuous at some critical value of y (see Fig. 
2). 

Similar modifications of the hysteresis loop are ob- 
served when changing T at constant y and p m Jh4|: at 
low T there are macroscopic avalanches on draining and 
on filling that disappear at two different critical temper- 
atures, Tf s {R) and T c des (i?); for higher T, both the ad- 
sorption and the desorption branches are continuous, and 
at Th(R) the two branches merge and the isotherm be- 
come fully reversible. This leads to the phase diagram 
illustrated in Fig. 1. Similar changes in the shape of the 
loop also occur when varying the porosity at constant 
yfL4|. This may account for the experimental observa- 
tions on aerogels [Q. 

Finally, we have also studied the connection between 
the above results and the characteristics of the free- 
energy landscape. Below the hysteresis line Th{R) (see 
Fig. 1), MFT predicts the appearance of a large num- 
ber of metastable states (minima of the grand-potential 
functional) in which the system can get trapped. Un- 
trapping only occurs under the influence of the exter- 
nal driving force, i.e., the change in /i. We again 
emphasize that this picture is valid provided the ther- 
mally activated processes for relaxation (processes that 
are for instance responsible for the creep motion of 
pinned interfaces [ft3|) are slower than the experimental 
timescale. The metastable states are responsible for the 
hysteresis loop discussed above as well as for the scanning 



curves and the inner subloops. Fig. 5 displays some of 
these states obtained by performing desorption scanning 
trajectories in the system without surface for y = 1.2 
(additional metastable states could be obtained by the 
procedure explained in Ref.Q). It is remarkable that 
the hysteresis loop obtained in the presence of the exter- 
nal surface appears as the true envelope of the scanning 
curves (compare with Fig. 2b), an observation that holds 
for all the cases studied. 

In summary, by combining mean-field density func- 
tional theory and a coarse-grained lattice-gas model, we 
have provided a comprehensive picture of hysteresis in 
capillary condensation of gases in disordered mesoporous 
solids. We have shown that the change of morphology of 
the hystereresis loop is related to the occurence of out- 
of-equilibrium phase transitions whose nature is different 
on adsorption and on desorption: this is due to the pres- 
ence of an interface between the porous matrix and the 
gas reservoir. 

The Laboratoire de Physique Theorique des Liquides 
is the UMR 7600 of the CNRS. 
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interface before pinning occurs diverges as a power law 
when [i — ► fj, c - 



